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QCD phase diagram for small densities from simulations at imaginary // 
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We present results on the QCD phase diagram for smaU densities without reweighting. Our simulations are 
performed with an imaginary chemical potential /ij for which the fermion determinant is positive. On an 8"^ x 4 
lattice with 2 flavors of staggered quarks, we map out the pseudo-critical temperature Tc{fJ,i). For /x/ /T < n/3, this 
is an analytic function whose Taylor expansion converges rapidly, with truncation errors smaller than statistical 
ones. The result is analytically continued to give the location of the pseudo-critical line for real /is S;500 MeV. 
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1. INTRODUCTION 

In view of heavy ion collision experiments a 
pressing task for lattice QCD is to map out the 
deconfinement transition in the {fi, r)-plane. The 
complexity of the fermion determinant at finite 
baryon density renders standard Monte Carlo 
techniques impossible [|l|. One way of circum- 
venting this 'sign problem' are multi-dimensional 
reweighting techniques, and a phase diagram has 
been presented last year [||. Reweighting gener- 
ally breaks down for large volumes and/or densi- 
ties, and one difficulty of the approach is to know 
when that happens. In both reweighting fac- 
tor and observables were expanded in a Taylor 
series in fi/T, and the first coefficients were mea- 
sured HJ. This method allows to simulate larger 
volumes, but a priori nothing is known about the 
error introduced by truncating the series. 

Here we present an alternative method avoiding 
reweighting altogether. We simulate at imaginary 
fj, for which there is no sign problem and arbitrar- 
ily large volumes are feasible. Observables are 
fitted by truncated Taylor series in /i/T, keeping 
full control over the associated systematic error. 
The series may then be continued to real /U, pro- 
vided the observable is analytic Q . This strategy 
has been successfully tested for screening masses 
|5j, here we extend it to the critical line itself. A 
detailed account of this work is given in [p| . 
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Figure 1. Phase diagram in the (/i/, T)-plane. 
The vertical line marks a Z{3) transition crossing 
the deconfinement transition. For T > the 
former is of first order. 



2. QCD AT IMAGINARY ^ 

Finite density QCD is symmetric under reflec- 
tion of fj,, Z{^) = Z{—^), for both real and imag- 
inary ^. For imaginary fi — ifii it is moreover 
periodic, Z{pi) = Z{pi + 27r/3) §. Hence, cer- 
tain shifts in fii are equivalent to Z(3) transfor- 
mations to different vacua. Consequently, Z(3) 
transitions occur at the critical values (/i/T)^ — 
(fc + l/2)27r/3. The Z(3) sectors can be identified 
by monitoring the phase of the Polyakov loop. 
Ref. M predicts these transitions to be of first or- 
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Figure 2. Probability distribution of the phase 
of the Polyakov loop for the critical value a/ij = 
7r/12. 

der for T > Tc and crossover for T < Tc, with the 
critical deconfinement temperature Tc- Schemat- 
ically, the phase diagram is shown in Fig. which 
is periodically repeated for larger /i/. 

Our simulations support this picture. Fig. ^ 
shows the phase of the Polyakov loop at {ji/TY 
for various lattice couplings/temperatures. A 
smooth distribution is found for low tempera- 
tures, whereas at high temperatures the two-state 
distribution signals a first order phase transition. 
Similar findings are reported for four flavors in 
The existence of a Z(3) transition with its 
endpoint joining on the deconfinement transition 
implies that the critical line of the latter be- 
comes non-analytic at that point, and hence only 
M/ ^ ""/S can be continued to real /i. 

Note that the deconfinement line as well as all 
observables even under refiections /i —fi are 
symmetric around the Z(3) transition. This is 
a feature that is not reproduced by reweighting 
from a Monte Carlo ensemble generated at /i = 0, 
since such an ensemble does not sufficiently probe 
the other Z(3) sectors §,||. 

3. CONTINUATION OF THE DECON- 
FINEMENT LINE 

One method to locate a phase transition is to 
look for the peak Xmax = xif^-c, Pc) of the suscep- 
tibilities X = VNt {{O - {0)f), where we have 
used the plaquette, the chiral condensate and the 
Polyakov loop as operators 0{x). Since simula- 
tions are always performed on finite volumes, sus- 
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Figure 3. Susceptibilities of the modulus of the 
Polyakov loop for various /i/ . 



ceptibilities are analytic functions over the whole 
parameter space of the theory. Non-analyticities 
associated with phase transitions develop only in 
the thermodynamic limit. 

For a given finite volume, the pseudo-critical 
line is then found by locating Xmax for every /i, 
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This is an implicit definition of the critical cou- 
pling Pdlj)- For analytic xil^^P) the implicit 
function theorem implies that also /3c (m) is ana- 
lytic for all fi. Moreover, the symmetry consider- 
ations of the last section ensure /3c(m) — /3c(— m), 
so that we can expand in /i/T = aNtfi to obtain 
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As long as, for imaginary = i/x/, the non- 
perturbative data are well described by a few co- 
efficients in this series, analytic continuation to 
real /i simply proceeds by the operation /i/ 
This provides us with the location of the pseudo- 
critical line on a finite volume. 

If there is a first order deconfinement transition 
ending in a critical point, then, as the volume 
is increased, /3c (m) approaches the line of phase 
transitions, while remaining pseudo-critical in a 
crossover regime. In principle, the nature of the 
line can be determined from the way the thermo- 
dynamic limit is approached, [(3c{V) — (3c{oo)) ~ 
V^"' , with o- = lorO<(7<l for a first or sec- 
ond order transition, respectively, and a = for 
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Figure 4. Location of the pseudo-critical decon- 
finement line as determined from plaquette sus- 
ceptibilities. The lines show fits of the form Eq. |^ 
to quadratic and quartic order. 
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Figure 5. Location of the deconfinement transi- 
tion. The error bar gives the uncertainty in Tc(0) 
used to set the scale, the dotted lines reflect the 
error on ci, Eq. ^j. 



a crossover. Our first simulations are on a single 
volume and this question is postponed to future 
work. 

Fig.|| shows a sequence of susceptibilities for 
the modulus of the Polyakov loop, with similar 
pictures for the plaquette and the chiral conden- 
sate. All three observables give statistically com- 
patible values /3c (a*/) • Note that the peak is mov- 
ing to larger Pc for growing /x/, as one would ex- 
pect for imaginary fj,, cf. Eq. || and Fig. [|. 

Fig. ^ shows Pc for different values of /i/. The 
data perfectly reproduce the expected symmetry 
of the partition function about the first Z(3) tran- 
sition point. The lines in the figure show fits of a 
quadratic and a quartic polynomial to the data, 
with x^/dof «0.6 and 0.7, respectively. The co- 
efficient of the quartic term is found to be zero 
within errors, and hence the two fits are statis- 
tically compatible: up to /x//T = 7r/3, and for 
the quark flavors and masses under consideration, 
the critical line is well described by the quadratic 
term only. 

After analytic continuation the result is trans- 
formed into physical units by means of the pertur- 
bative two-loop expression for the scale parameter 
qAl and using Tc{n — 0) from the literature 
to set the physical scale. We thus arrive at the 
phase diagram of Fig. with the pseudo-critical 
line given by (/xs = 3/z) 



TcjfJ-B) 
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0.00563(38) 



This result is consistent with those obtained by 
the reweighting approaches [^,^ , for which it pro- 
vides a crucial, controlled check. An important 
task ahead is to perform a finite size scaling anal- 
ysis of the continued critical line, in order to de- 
termine the order of the transition as well as to 
locate a possible critical point. 
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